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The first mitochondrial genome for the butterfly family Riodinidae 
(Abisara fylloides) and its systematic implications 
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Abstract: The Riodinidae is one of the lepidopteran butterfly families. This study describes the complete mitochondrial genome of 
the butterfly species Abisara fylloides, the first mitochondrial genome of the Riodinidae family. The results show that the entire 
mitochondrial genome of A. fylloides is 15301 bp in length, and contains 13 protein-coding genes, 2 ribosomal RNA genes, 22 
transfer RNA genes and a 423 bp A+T-rich region. The gene content, orientation and order are identical to the majority of other 
lepidopteran insects. Phylogenetic reconstruction was conducted using the concatenated 13 protein-coding gene (PCG) sequences of 
19 available butterfly species covering all the five butterfly families (Papilionidae, Nymphalidae, Peridae, Lycaenidae and 
Riodinidae). Both maximum likelihood and Bayesian inference analyses highly supported the monophyly of Lycaenidae+Riodinidae, 
which was standing as the sister of Nymphalidae. In addition, we propose that the riodinids be categorized into the family Lycaenidae 


as a subfamilial taxon. 


Keywords: Abisara fylloides; Mitochondrial genome; Riodinidae; Systematic implication 


The typical metazoan mitochondrial genome 
(mitogenome) contains 37 genes, including 13 protein- 
coding genes (PCGs), 2 rRNA genes and 22 tRNA genes, 
and a non-coding area (i.e., the control region or the 
A+T-rich region) (Wolstenholme, 1992; Boore, 1999). 
Maternal inheritance, lack of recombination and an 
accelerated evolutionary rate compared with the nuclear 
genome have all contributed to the increased use of 
mitogenomes, which is one of the key methods in fields 
such as phylogenetics, comparative and evolutionary 
genomics, molecular evolution and population genetics 
(Ballard & Whitlock, 2004; Simonsen et al, 2006). At 
present, mitochondrial genomes have already been 
determined in a variety of insect groups covering nearly 
200 species. However, reported complete mitogenomes 
are relatively scarce for lepidopterans and especially for 
butterflies. To our knowledge, as of October 2012 only 
about 20 butterfly species covering five butterfly families 
(Table 2) have been reported or deposited into the 
GenBank, but only one butterfly family, the Riodinidae, 
still lack corresponding data. 

The phylogenetic position and taxonomic ranking of 
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the butterfly family Riodinidae among butterfly lineages 
are still controversial issues among entomologists. Some 
scholars suggest that the riodinids are closely related to the 
lycaenids, considering the similarities in morpholo gical 
character, behavior, and host plants between the two (slug- 
like larvae, pupa contigua, ants associated) (Ackery, 1984; 
Chou, 1998; de Jong et al, 1996; Ehrlich, 1958; Scott, 
1985). Moreover, they are usually classified into the 
Lycaenidae family as a subfamilial taxon. Some consider 
the riodinids a unique family parallel to the Lycaenidae 
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family (Campbell et al, 2000; de Jong et al, 1996; 
Kristensen, 1976; Shou et al, 2006), and others propose 
that the riodinids are more closely related to the 
nymphalids than to the other butterfly groups, such as the 
lycaenids (Martin & Pashley, 1992; Robbins, 1987, 1988). 
This study sequenced the first complete 
mitochondrial genome of the Abisara fylloides, a 
representative species of the family Riodinidae, by long 
PCR and primer walking techniques. Its genetic structure 
was preliminarily compared with those of other available 
butterfly species. The maximum likelihood (ML) and 
Bayesian inference (BI) phylogenetic trees of the A. 
fylloides and other available butterfly representative 
species were reconstructed based on concatenated DNA 
sequences of the 13 protein-coding genes (PCGs), and aim 
to clarify their phylogenetic relationships and further 
provide new information about the structure, organization 
and molecular evolution of the lepidopteran mitogenomes. 


MATERIALS AND METHODS 


Sample collection and DNA extraction 
Adult individuals of A. fylloides were collected in 


Jinghong, Yunnan Province, China in August 2006. After 
collection, sample specimens were preserved in 100% 
ethanol immediately and stored at —20 °C before DNA 
extraction (specimen No. ZHF07). Total Genomic DNA 
was isolated using the proteinase K-SiO, method as 
described by Hao et al (2005). 
PCR amplification and sequence determining 

The multiple sequence alignments were conducted 
using the software Clustal X 1.8 based on the mitogenome 
sequences of Coreana raphaelis, Artogeia melete, Troides 
aeacus available from GenBank and those of Argyreus 
hyperbius, Acraea issoria, Calinaga davidis, Pieris rapae 
determined in our laboratory (Thompson et al, 1997). The 
long PCR primers, which may cover the whole 
mitogenome, were designed according to the conserved 
regions by the software Primer premier 5.0 (Singh et al, 
1998) (Table 1). Seven short fragment sequences (500-700 
bp) of сохі, cox2, cox3, cytb, пайї, rrnL and rmS were 
amplified using insect universal primers (Caterino & 
Sperling, 1999; Simmons & Weller, 2001; Simon et al, 
1994). All the primers were synthesized by the Shanghai 
Sheng gong Biotechnology Co. Ltd. 


Table 1 List of PCR primers used in this study 





Genes Forward primers (5) Reverse primers (3’) ы C) 
сох1* GGTCAACAAATCATAAAGATATTG TAAACTTCAGGGTGACCAAAAAT 50.0 
сох1-сох2 ТТАТТТСТАТСАССССТАС ATAGCAGG AAGATTGTTC 47.5 
cox2* GAGACCATTACTTGCTTTCAGTCACT CTAATATGGCAGATTATATGTATGG 49.5 
cox2-cox3 TTTTATTGCTCTTCCATCT TTATTCCTCATCGTAATCC 48.5 
cox3* TATTTCAATGATGACGAGAT CAAATCCAAAATGGTGAGT 49.8 
cox3-nad5 TTTATAGCAACAGGATTTC CATCAACTGGTTTAACTTT 45.5 
nad5 AAAACTTCCAGAAAATAATCTC TTGCTTTATCTACTTTAAGACA 46.5 
nad5-cytb AATTATACCAGCACATAT TTATCGACTGCAAATC 47.1 
cytb* TATGTACTACCATGAGGACAAATAT ATTACACCTCCTAATTTATTAGGAAT 47.0 
cytb-nadl TCCTGCTAACCCTTTAGTCA AGGTAGATTACGGGCTGTT 48.0 
nad1* CGTAAAGTCCTAGGTTATATTCAGATCG ATCAAAAGGAGCTCGATTAGTTTC 52.0 
nadl-rrnL AGCCCGTAATCTACCTAA TAAGACGAGAAGACCCTAT 47.0 
rmL* CGCCTGTTATCAAAAACAT CCGGTCTGAACTCAGAT 45.5 
rrnL-rrnS AGACTATTGATTATGCTACCT TAAGAATCTAATGGATTACAA 46.5 
rm s* CTTCTACTTTGTTACGACTTA T AATTTTGTGCCAGCAGTTG 50.0 
rrnS-nad2 AGAGGGTATCTAATCCGAGTTT TGGCTGAGAATTAAGCGATA 49.5 
nad2-coxl ATACAGAAGCAGCATTA AGAAGGAGGAAGAAGTCAA 52.0 


*: universal primer. 


Seven partial gene sequences were initially 
sequenced under the following conditions: an initial 
denaturation at 94 °С for 5 minutes, then denaturation at 
94 °С for 1 minute for a total of 35 cycles; annealing at 
45—55 °С for 1 minute and extension at 72 °С for 2 
minutes plus 30 seconds; final extension at 72 °С for 10 
minutes. Long PCRs were performed using TaKaRa LA 
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Taq polymerase with the following cycling parameters: 
an initial denaturation for 5 minutes at 95 ?C; followed 
by 30 cycles at 95 °С for 55 seconds, 45—55 °С for 2 
minutes, 68 ?C for 2 min and 30 seconds; and a 
subsequent final extension step of 68 °С for 10 minutes. 
The PCR products were separated by electrophoresis 
in a 1.296 agarose gel and purified using the DNA gel 
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extraction kit (TaKaRa). All PCR fragments were 
sequenced directly after purification with the QIA quick 
PCR Purification Kit reagents (QIAGEN). Internal primers 
were applied to complete sequences by primer walking 
(detailed primer information will be provided upon request). 
All fragments were sequenced for both strands. 


Data analysis 

We used DNASIS MAX (Hitachi) for sequence 
assembly and annotation. Protein-coding genes and 
rRNA genes were identified by sequence comparison 
with other available insect mitochondrial sequences. The 
tRNAs were identified by tRNAscan-SE v.1.21 (Lowe & 
Eddy, 1997). The putative tRNAs, which were not found 
by tRNAscan-SE, were identified by a sequence 
comparison of A. fylloides with the other lepidopteran 
tRNAs. PCGs were aligned with the other available 
lepidopteran mitogenomes using DAMBE software (Xia 
& Xie, 2001). The tandem repeats in the A+T-rich region 
were predicted using the Tandem Repeats Finder online 
(http://tandem.bu.edu/trf/trf.html) ^ (Benson, 1999). 
Nucleotide composition was calculated using PAUP 
4.0b10 (Swofford, 2002). The mitogenome sequence 
data have been deposited in GenBank under the 
accession number HQ259069. 


Phylogenetic analysis 

Phylogenetic analyses were performed on 19 
representative species including A. fylloides, covering all 
the six families of butterflies. The multiple aligning of 
the concatenated nucleotide sequences of the 13 
mitochondrial PCGs of the 19 species (Table 2) was 
conducted using ClustalX 1.8. The phylogenetic trees 
were reconstructed with the maximum likelihood (ML) 
and Bayesian inference (BI) methods, using the moth 
species Adoxophyes honmai (GenBank accession 
number of mitogenome: NC014295) as the outgroup. In 
both phylogenetic analyses, the third codon position of 
all the sequences was excluded. The ML analyses were 
conducted in PAUP 4.0b10 by using TBR branch 
swapping (10 random addition sequences) as a search 
method. The model GTR+I+T was selected as the best fit 
model using Modeltest 3.06 (Posada & Crandall, 1998) 
under the AIC scores, and the bootstrap values of the ML 
tree were evaluated via the bootstrap test with 1 000 
iterations. The Bayesian analysis was performed using 
MrBayes 3.1.2 (Huelsenbeck & Ronquist, 2001) with the 
partitioned strategy (13 partitions: cox], cox2, cox3, atp8, 
atp6, nad1, nad2, nad4, nad4L, nad5, nad6 and cytb), and 
the best substitution model for each partition was 
selected as in the ML analysis. The MCMC analyses 
(with random starting trees) were run with one cold and 
three heated chains simultaneously for 1,000,000 
generations sampled every 100 generations with a burn- 
in of 25% until the average standard deviation of split 
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frequencies to be less than 0.01, which means that 
convergence was reached. 


RESULTS AND DISCUSSION 


Genome structure and organization 

The complete mitogenome is 15 301 bp in length, 
encodes 37 genes in all. It contains 13 protein, 22 tRNA, 
2 rRNA genes and a non-coding high A+T content region 
(Figure 1, Table 3). Its structure and organization are 
identical to those of the majority of other lepidopterans 
(Bae et al, 2004; Cha et al, 2007; Cameron & Whiting, 
2008; Hao et al, 2012; Hong et al, 2008, 2009; Hu et al, 
2010; Kim et al, 2009b; Ji et al, 2012; Junqueira et al, 
2004; Wang et al, 2011), though a few lepidopterans, 
such as three Thitarodes species, were reported to 
possess the ancestral gene order trnI-trnQ-trnM instead 
of the trnM-trnI-trnQ (Cao et al, 2012) . 

Eight overlapping sequences totaling 61 bp are 
located throughout the A. fylloides mitogenome, with 
size ranging from 2 to 35 bp, of which the longest (35 bp) 
is located between the cox2 and the tRNA genes. In 
addition, 17 intergenic spacers ranging from 1 to 45 bp in 
length are found in the mitogenome. Among these 
spacers, the longest is located between the КМА" and 
nad2 genes, the other 16 spacers are scattered throughout 
the whole genome (Figure 1, Table 3). Most of these 
spacer regions are arranged relatively compactly 


Abisara fyuoides 


mtDNA 15301 bp 





Figure 1 Circular map of the mitochondrial genome of 
Abisara fylloides 

Gene names not underlined indicate the direction of transcription from left to 

right and those underlined indicate right to left. Transfer RNA genes encoded 

by H and L strands are shown outside and inside the circular gene map, 

respectively. Transfer RNA genes are indicated by the IUPAC-IUB single letter 

amino acid codes, while L1, L2, S1, S2 represent tRNA-Leu(UUR), tRNA- 


Leu(CUN), tRNA-Ser(AGN) and tRNA-Ser(UCN), respectively. 
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Table 3 Organization of the Abisara fylloides mitochondrial genome 





Gene Direction Position Size (bp) Intergenic length” Start codon Stop codon 
tRNAM* F 1-67 67 0 
tRNA! F 68-131 64 -3 
tRNA?" R 129-197 69 45 
nad2 F 243-1 256 1014 1 ATT TAA 
tRNA"? F 1 258-1 324 67 -8 
tRNA R 1317-1382 66 3 
tRNA’ R 1386-1 451 66 6 
сох1 Е 1 458-2 989 1 532 3 CGA T 
tRNALUOLR) F 2 993-3 059 67 0 
cox2 F 3 060-3 738 679 3 ATG T 
tRNA F 3 736-3 806 71 
tRNA‘? F 3 807-3 871 65 
atp8 F 3 872-4 033 162 -7 ATC TAA 
atp6 F 4 027-4 714 688 -2 ATG TAA 
сох3 Е 4713-5 501 789 2 ATG TAA 
tRNAGY F 5 504—5 569 66 0 
nad3 F 5 570—5 923 354 3 ATT TAA 
tRNAA" F 5 927—5 997 71 3 
tRNAA® F 6 001-6 064 64 1 
tRNA n F 6 066—6 131 66 16 
ЕТА АЧА» Е 6 148—6 208 61 1 
tRNAS" F 6 210-6 275 66 -2 
tRNA’ R 6 274-6 340 67 0 
nad5 R 6 341-8 108 1 768 ATT T 
tRNA" R 8 109-8 179 71 
nad4 R 8186-9 518 1333 E ATG T 
nad4L R 9 517—9 805 289 2 ATG TAA 
tRNA!” F 9 808—9 870 63 0 
tRNA"? R 9 871-9 937 67 5 
nad6 F 9 943-10 467 525 3 ATA TAA 
cytb F 10 471-11 622 1152 -2 ATG TAA 
tRNASEHUCN) F 11 621-11 685 65 17 
пайї R 11 703-12 641 939 1 ATG TAA 
tRNAP CUN R 12 643-12 710 68 0 
rrnL R 12 711-14 044 1334 0 
tRNA I R 14 045-14 107 63 0 
rrnS R 14 108-14 878 771 0 
A+T-rich region 14 879-15 301 423 


`: In the column intergenic length, the positive number indicates interval base pairs between genes, while the negative number indicates the overlapping base 


pairs between genes. 


compared with other insect mitogenomes (Cameron & 
Whiting, 2008; Hao et al, 2012; Hong et al, 2009; Hu et 
al, 2010; Kim et al, 2009b; Ji et al, 2012; Junqueira et al, 
2004; Wang et al, 2011). 

The A. fylloides A+T-rich region is flanked on one 
side by the rmS and on the other side by the tRNA“ 
genes. This region exhibits the highest A+T content 
(91.0%) (Table 4), and spans 423 bp (Table 3, Table 4). A 
sequence analysis of the A+T-rich region revealed that it 
contained some structures typical of other lepidopteran 
mitogenomes: (1) at 20 bp downstream of the small 
subunit rRNA gene, there is a structure including a 
motif ‘АТАСА? which is very well conserved in all 
sequenced lepidopteran insects, and a 18-bp polyT 
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stretch, both of which have been suggested as the origin 
of minority or light strand replication (On) and to play a 
regulatory role (Kim et al, 2009a; Lutz-Bonengel et al, 
2004; Saitou et al, 2005; Yukuhiro et al, 2002). (2) 
Between the sites 15 151 and 15 168 there is a 
microsatellite-like (TA) element, which is also found 
in the majority of other lepidopterans; (3) There is 
another motif “ATTTA” of unknown function located 
from 15 139 to 15 143 upstream of ће (TA)s, which is 
also typical of the other lepidopterans. In addition, to 
our great surprise, an unexpected short microsatellite- 
like repeating region (ТА), downstream of the (ТА) 
was detected in the А+Т-гїсһ region, and this has not 
been reported in any other lepidopterans. 
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Table 4 Nucleotide composition and skewness in different regions of the Abisara fylloides mitogenome 


Nucleotide composition (%) 








Region в AT-skew GC-skew 
(bp) T б А G А+Т 

Whole genome 15 301 41.7 11.3 39.5 7.5 812 —0.027 —0.202 
Major-strand PCGs 6 927 44.5 12.1 34.1 9.3 78.6 —0.132 —0.131 
Minor-strand PCGs 4329 47.6 6.3 34.1 12.0 81.7 —0.165 0.311 
Whole PCGs 11 224 45.8 9.9 34.0 10.3 79.8 —0.148 0.020 
Whole tRNA 1 461 40.0 TA 41.5 10.7 81.5 0.018 0.163 
rmL 1 334 38.4 48 47.0 9.8 85.4 0.101 0.342 

rms 771 39.9 4.9 45.7 9.5 85.6 0.068 0.319 
A+T-rich region 423 48.9 4.5 42.1 4.5 91.0 —0.075 0.000 


Base composition bias 

The base composition of the А. fylloides 
mitogenome shows an A and T bias (Table 4). The 
whole А+Т content of the mitogenome 15 up to 81.295, 
ranging from that of Argyreus hyperbius (80.81%) to 
that of Coreana raphaelis (82.66%). Like most other 
metazoan mitogenomes, the A+T content of the A+T- 
rich region, which is located between the rrnS and 
tRNA М“ genes, is the highest (91.0%) in all known 
butterfly species except for P. rapae (Pieridae) to date. 
The base contents of A, T, C, G are 39.5%, 41.7%, 
11.3%, 7.5%, respectively, indicating a relatively 
higher A+T content (81.2%). These phenomena 
commonly exist in the protein-coding genes, which 
have a relatively lower A+T content (79.8%), and the 
tRNA and rRNA genes in insects. 


Protein-coding genes 

The sequences of the 13 Abisara fylloides PCGs 
are 11 224 bp in length, including 3 730 codons 
(excluding termination codons). Twelve of the 13 
PCGs use standard ATN as their start codon except for 
the coxl gene, and eight of these 12 PCGs begin with 
ATA or ATG (Methionine) and the other four begin 
with АТТ or ATC (Isoleucine) (Table 3). The start 
codons for coxl gene of lepidopteran insects have 
usually been a controversial issue. In general, there is 
no typical start codon, especially for the сох1 gene, 
which usually uses CAG (R) as start codon in insects, 
including the lepidopterans. However, some scholars 
reported unusual start codons, such as the 
trinucleotide TTG (Bae et al, 2004; Hong et al, 2008), 
ACG (Lutz-Bonengel et al, 2004), GCG (Nardi et al, 
2003), the tetranucleotide ATAA, ATCA and ATTA 
(Clary & Wolstenholme, 1983; de Bruijn,1983; Kim et 
al, 2006), and the hexanucleotides TATTAG (Flook et 
al, 1995), TTTTAG (Yukuhiro et al, 2002), TATCTA 
(Coates et al, 2005), ATTTAA (Beard et al, 1993; 
Mitchell et al, 1993) for coxl in some other insect 
species. However, the CGA is present as a conserved 
region for all lepidopteran insects reported, and thus 
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we tend to consider that CGA is the сох1 start codon 
for A. fylloides as Kim et al (2009b) suggested. As for 
stop codons, 9 of the 13 PCGs use standard TAA 
except for the coxl, cox2, nad4 and nad5 genes, all of 
which terminate at a single thymine (Table 3), and this 
case was also found in all other lepidopterans reported 
to date. For more details on this phenomenon, please 
refer to the discussions of Kim et al. (2010). 

The PCG amino acid sequence variation analysis 
showed that there were 3 755 homologous sites, of 
which 1 964 are conserved, 1 791 are variable, and 1 
128 are parsimony informative. Among the twenty 
amino acids in the 13 PCGs, six (Leu, Met, Ile, Phe, 
Asn, and Tyr) were used more frequently than the 
others, and their usage frequencies were higher than 
the average, whereas the other 14 amino acids less 
used (Table 5). 


rRNA and tRNA genes 

The large subunit rRNA and small subunit rRNA 
genes of the A. fylloides are 1 334 and 771 bp in length, 
respectively. As in other lepidopterans, these two genes 
are located between tRNA!” and tRNA, and 
between tRNA I and A+T-rich region respectively 
(Cameron & Whiting, 2008; Hao et al, 2012; Hu et al, 
2010; Kim et al, 2009a; Salvato et al, 2008; Wang et al, 
2011; Yang et al, 2009). 

The A. fylloides mitogenome harbors 22 tRNA 
genes, which are scattered throughout the whole genome 
and ranged in length from 61 to 71 bp. Except for the 
tRNA SGN, which lacks the DHU loop, all tRNAs are 
shown to be folded into the cloverleaf secondary 
structures, within which all amino acid acceptor stems 
have 7 base pairs, and all anticodon stems have 5 base 
pairs. This was also found in all the other lepidopterans 
determined to date (Cameron & Whiting, 2008; Hao et al, 
2012; Hu et al, 2010; Kim et al, 2009a). 

A total of 27 pairs of base mismatches were 
detected in all the predicted tRNA secondary structures, 
among which 17 are GU, 6 are UU, 2 are AA, and 2 are 
AC. The AA mismatches occur at the anticodon stem of 
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Table 5 Codon usage of the protein-coding genes of the Abisara fylloides mitogenomes 








Codon N (RSCU) Codon N (RSCU) Codon N (RSCU) Codon N (RSCU) 
UUU (F) 356 (1.86) UCU (S) 108 (2.59) UAU (Y) 175 (1.83) UGU (С) 28 (1.81) 
UUC (F) 26 (0.14) UCC (S) 16 (0.38) UAC (Y) 16 (0.17) UGC (C) 3 (0.19) 
UUA (L) 488 (5.30) UCA (S) 94 (2.25) UAA (*) 0 (0) UGA (W) 84 (1.87) 
UUG (L) 12 (0.13) UCG (S) 2 (0.05) UAG (*) 0 (0) UGG (W) 6 (0.13) 
CUU (L) 34 (0.37) CCU (P) 70 (2.33) CAU (H) 62 (1.85) CGU (R) 14 (1.04) 
CUC (L) 3 (0.03) CCC (P) 9 (0.3) CAC (H) 5 (0.15) СОС (В) 0 (0) 
CUA (L) 14 (0.15) CCA (P) 38 (1.27) CAA (О) 57 (1.93) CGA (R) 37 (2.74) 
CUG (L) 1 (0.01) CCG (P) 3 (0.1) CAG (Q) 2 (0.07) CGG (R) 3 (0.22) 
AUU (D 424 (1.88) ACU (T) 74 (2.1) AAU (N) 251 (1.81) AGU (S) 25 (0.6) 
AUC (D 27 (0.12) ACC (T) 7 (0.2) AAC (N) 26 (0.19) AGC (S) 0 (0) 
AUA (M) 283 (1.81) ACA (T) 59 (1.67) AAA (К) 97 (1.92) AGA (S) 89 (2.13) 
AUG (M) 29 (0.19) ACG (T) 1 (0.03) AAG (K) 4 (0.08) AGG (S) 0 (0) 
GUU (V) 74 (2.39) GCU (A) 73 (2.52) GAU (D) 55 (1.83) GGU (G) 56 (1.17) 
GUC (V) 4 (0.13) GCC (A) 8 (0.28) GAC (D) 5 (0.17) GGC (G) 1 (0.02) 
GUA (V) 45 (1.45) GCA (A) 33 (1.14) GAA (E) 65 (1.78) GGA (G) 100 (2.08) 
GUG (V) 1 (0.03) GCG (A) 2 (0.07) GAG (E) 8 (0.22) GGG (G) 35 (0.73) 





n: frequency of codon used; RSCU: relative synonymous codon usage; *:stop codon. 


Stop codons were excluded in total codon counts. 


tRNA'^ and the TyC loop of tRNA"; the AC 
mismatches occur at the amino acid acceptor stem of 
tRNA " and the anti codon stem of tRNA"), Among 
the 6 UU mismatches, two occur at the amino acid 
acceptor stems of tRNA г and tRNA F, two at the 
anti codon stems of {ВМА (ЧУ and ВМА", and 
another two at the anti codon stem of tRNA?* CV, 
respectively (Figure 2). 


Phylogenetic analysis 

At present, for the phylogenetic positions of the 
riodinids within —papilionid butterflies, most 
morphological studies place them as most closely related 
to the lycaenids and identify the nymphalids as the 
closest relatives to this riodinid+lycaenid clade (de Jong 
et al, 1996; Ehrlich & Ehrlich, 1967; Kristensen, 1976; 
Scott & Wright, 1990). These relationships have been 
inferred using a variety of phylogenetic methods and are 
supported by a number of adult, larval and pupal 
synapomorphies. Additionally, molecular (DNA 
sequence of the mitochondrial МАРНІ gene) or 
molecular plus morphological evidence also result in a 
monophyletic interpretation of the Riodinidaet+ 
Lycaenidae, and their sister relationship to the 
Nymphalidae (Wahlberg et al, 2005; Weller et al, 1996). 
However, based on a cladistic analysis of four foreleg 
characters with nine character states, Robbins (1988) 
suggested that the Riodinidae are more closely related to 
the Nymphalidae than to the Lycaenidae, and this result 
is supported by the nuclear 28S rRNA gene sequence 
data (Martin & Pashley, 1992). 

There are two opinions regarding the taxonomic 
rank of riodinids. First, some previous studies suggest 
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that the riodinids should be categorized into the family 
lycaenids as a subfamilial taxon in light of their 
morphological characters (Chou, 1998; de Jong et al, 
1996; Ehrlich, 1958; Kristensen, 1976; Scott & Wright, 
1990), and this opinion is supported by the molecular 
studies of Zou et al (2009) and Hao et al (2007). Other 
studies postulate that the riodinids should be classified as 
a separate family parallel to Lycaenidae (Harvey, 1987; 
Martin & Pashley, 1992; Robbins, 1988; Weller et al, 
1996), and this view is supported by molecular 
phylogenetic studies based on data from the wingless 
gene by Campbell et al (2000) and the combined analysis 
of sequences of the nuclear Ef-la, wingless, and 
mitochondrial COI genes by Wahlberg et al (2005). 

The ML and Bayesian trees of this study (Figure 3) 
showed that all the butterfly taxa in this study did not 
form a monophyletic unit, and a similar case was 
reported by Hao et al (2012). Nonetheless, both the ML 
and BI trees indicated that all the butterfly species were 
grouped into five distinct lineages: 1) the Papilionidae, 
including papilionids and parnassids; 2) Hesperiidae; 3) 
Pieridae; 4) Nymphalidae; 5) Lycaenidae+Riodinidae. 
The monophyly of Lycanidae + Riodinidae was strongly 
supported with a 100% bootstrap value in ML, and with 
1.00 posterior probability value in BI. Thus, considering 
the results of Heikkilà et al (2012), which indicate the 
monophylies of lycaenids and riodinids, it is reasonable 
to propose that the two groups may be sisters, though the 
taxa sampling of riodinids in this analysis is extremely 
limited. Additionally, based their congruent genetic 
divergences compared with those between other butterfly 
subfamilies, the riodinids should be categorized into the 
Lycaenidae family as a subfamilial taxon. 
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Figure2  Predicated clover-leaf secondary structures for the mitochondrial tRNA genes of Abisara fylloides 
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Figure3 The Bayesian inference (BI) and maximum likelihood (ML) phylogenetic trees of main butterfly lineages based on 13 


protein-coding gene sequences (Numbers on each node correspond to the posterior probability values of the BI analysis and the ML 


bootstrap percentage values for 1 000 replicates of ML analysis) 
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Abstract: There is some discrepancy in the classification of different species of Hylopetes, particularly regarding systematic status of 
H. electilis and H. phayrei and their relationship to other species. In the present study, for the first time we have brought together six 
of the nine Hylopetes species and performed statistical analysis of 14 measurable cranial variables, analyzing in total 89 specimens, 
including H. electilis, H. alboniger, H. phayrei, H. lepidus, H. spadiceus, and H. nigripes. Both univariate and multivariate analysis 
results indicate that H. electilis can not only be obviously distinguished from H. phayrei, but also clearly differs from the other four 
Hylopetes species. These results sustain the contention that H. electilis is neither a synonym nor subspecies of H. phayrei, but should 
be considered a distinct and valid species. Subsequently, a straightforward discussion on the biogeography of Hylopetes in 


southeastern Asia gives further insight into the differentiation and variety of species belonging to this genus. 


Keywords: Hylopetes; Morphology; Cranial variables; Statistical analysis; Species 


The genus Hylopetes (Thomas, 1908) mainly 
inhabits mountain forests in Asia, ranging from the 
Himalayas to the Greater Sunda Islands, and including 
areas of Nepal, China, Indochina, Philippines, Malaysia, 
and Indonesia (Corbet & Hill, 1992; Hoffmann et al, 
1993; Nowak, 1999; Thorington & Hoffmann, 2005). 
Initially, Hylopetes was regarded as a subgenus of 
Sciuropterus (Thomas, 1908), but, based on Pocock’s 
(1923) studies on the morphology of baculum, it was 
elevated to a full genus. This view, that the arrow-tailed 
flying squirrel Hylopetes in the eastern extreme of the 
Himalayas and southeastern Asia is a valid genus from the 
small Kashmir flying squirrel Eoglaucomys (A. H. Howell, 
1915) in the western Himalayas, has been sustained by 
subsequent morphological and molecular studies (Nowak, 
1999; Oshida et al, 2004; Thorington et al, 1996). 

Based on several dental and cranial characteristics 
as well as external structures, different taxonomists have 
described varying numbers of Hylopetes species (Allen, 
1940; Corbet & Hill, 1992; Ellerman, 1940; Ellerman & 
Morrison-Scott, 1950; Hoffmann et al, 1993; Nowak, 
1999; Thorington & Hoffmann, 2005). To date, the 
maximal number of Hylopetes species was put at 13 
(Ellerman, 1940), though the latest accepted result is 9 


Science Press 


(Thorington & Hoffmann, 2005). 

Allen (1940) described two subspecies of Pteromys 
(Hylopetes) alboniger (Allen, 1925), and listed the Hainan 
flying squirrel as Pteromys electilis (Allen, 1925). 
Ellerman (1940) listed 13 Hylopetes species, including H. 
platyurus (Jentink, 1890), H. spadiceus (Blyth, 1847), H. 
phayrei (Blyth, 1859), H. alboniger (Hodgson, 1836), and 
H. nigripes (Thomas, 1893) etc., but the Hainan flying 
squirrel from Namfong, Hainan Island, China was 
regarded as Petinomys electilis (Allen, 1925). Ellerman & 
Morrison-Scott (1950) accepted the same results of 
Petinomys electilis as Ellerman (1940), however, they 
listed only five Hylopetes species, including H. spadiceus, 
H. phayrei, H. alboniger, H. fimbriatus (Gray, 1837), and 
Н. sagitta (Linnaeus, 1766). The last two species were 
emended in later studies: H. fimbriatus was listed as 
Eoglaucomys fimbriatus and H. sagitta was replaced by H. 
lepidus (Horsfield, 1822) in the results of Thorington & 
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Hoffmann (2005). Corbet & Hill (1992) regarded Н. р. 
electilis as subspecies of H. phayrei, and Sciuropterus 
platyurus was accepted as the synonym of H. lepidus 
(Horsfield, 1822), resulting in their listing of ten Hylopetes 
species, except for H. fimbriatus and H. baberi (Blyth, 
1847), which were arranged as Eoglaucomys fimbriatus, 
and H. platyurus (Thorington & Hoffmann, 2005). The 
rest of the eight Hylopetes species were the same as the 
results of Thorington & Hoffmann (2005). 

Despite lengthy and protracted debate, there 
remains several disagreements and uncertainties as to the 
taxonomy of this polymorphic genus. Among the most 
contentious include H. electilis in Hainan Island, China 
and H. alboniger, H. phayrei, H. platyurus, H. spadiceus, 
and H. nigripes, which have had a complicated taxon- 
omic history due to intraspecific geographic variation in 
their distributions across Asia (Corbet & Hill, 1992; 
Hoffmann et al, 1993; Rasmussen & Thorington, 2008; 
Roberts, 1997; Thorington & Hoffmann, 2005). Despite 
the controversies, to date no attempt has been made to 
investigate the differences in measurable skull variables 
among these groups. 

Morphometric data has been used to evaluate 
cranial, dental, and body measurements of various wild 
populations of mammals, and the resulting statistic 
analysis of variations on the morphological level has 
been proved useful in detecting patterns of geographic 
variations and delimiting intra- or  inter-specific 
evolutionary units (Li et al, 2008; Li et al, 2012; Munoz 
& Perpinan, 2010; Slabova & Frynta, 2007; Zelditch et al, 
2004). Because the specimens of Hylopetes are 
undoubtedly exiguous, we based our analysis on 89 
specimens of 6 Hylopetes species. Principal components 
analysis (PCA) was performed to compare skull 
measurable variables among them, and furthermore, 
special attention was paid to the taxonomic status of 
Hainan flying squirrel from Namfong, Hainan Island, 
China, which has been listed as either a synonym or 
subspecies of H. phayrei (Corbet & Hill, 1992; 
Hoffmann et al, 1993; Thorington & Hoffmann, 2005). 


MATERIALS AND METHODS 


Samples 

A total of 89 specimens of Hylopetes were examined 
and used for statistical analysis, including 39 females, 44 
males, and 6 specimens of unknown sex (Appendix A). 
Adults were classified based on phalange ossification and 
dental patterns. All specimens used in this study were 
deposited in the American Museum of Natural History 
(New York, AMNH); the National Museum of Natural 
History (Washington DC, USNM); the Institute of Zoology, 
Chinese Academy of Sciences (IOZ, CAS) (Beijing, China); 
and the Kunming Natural History Museum of Zoology, 
Kunming Institute of Zoology (KNHMZ, KIZ, CAS) 
(Kunming, China) In total, 14 measurable cranial 
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variables (Figure 1) taken with a digital caliper to the 
nearest 0.01 mm were used in statistic analysis. 











Figure 1 Measurements of 14 cranial variables used in the study 


Data analysis 

Statistical analysis of the morphometric study, which 
included sexual dimorphism, multiple comparisons, and 
principal component analysis (PCA) were performed 
using SPSS 11.0 (SPSS Inc., Chicago, USA). First, all 
variables were transformed into natural logarithms to 
eliminate the bias effect on large measurements in 
Statistical analysis (D’elia & Pardinas, 2004). Then, 
Kaiser-Meyer-Olkin Measure of Sampling Adequacy 
(KMO) and Bartlett’s Test were performed in order to 
test if the data were fit for analysis. Statistical differences 
were considered significant where P<0.05. 

The Tests of Equality of Group Means (by Wilks’ 
Lambda) was used to assess the sexual dimorphism of 
each group. In order to evaluate variations between 
samples, multiple comparisons (by Least-Significant 
Difference) between taxa were performed for all 14 cranial 
measurements. The associations between cranial meas- 
urable characters and species were assessed by multi- 
variate analysis of the PCA. PCA based on the variance- 
covariance matrix of the log-transformed variables was 
performed to identify variables that account for maximum 
variation in the data set, and to represent distances 
between major groups in order to assess the specific 
relationships among the individuals. The eigenvector 
scores describing the relative significance of each variable 
to principal components were used to compare the cranial 
morphological similarities апа differences. The PCA 
scatter-plot visually represents the variation among 
different individuals of the samples (Figure 2). 
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Figure 2 Plots of six Hylopetes species samples on the 
principal components factors 1 and 2 


RESULTS 


The means and standard deviation (SD) of 14 
measurable cranial variables of Hylopetes are prese- 
nted in Table 1. The results of КМО 15 0.946 and that 
of P of Bartlett's Test of Sphericity is 0.000, together 


confirming that the data garnered in the morphometric 
study were fit for PCA analysis. The Tests of Equality 
of Group Means for 83 samples showed no sexual 
dimorphism in all 14 cranial variables in the 6 Hylopetes 
species (Table 2). Table 3 presents the results of multiple 
comparisons on all 14 cranial variables between the 
species of Hylopetes. 

Multivariate analysis was based on principal 
components analysis. In PCA, eigenvalues for the first 
three principal components were 12.161, 0.740 and 
0.541 respectively, accounting for 96.014% of the total 
variance (Table 4). The first principal component factor 
(РСІ) accounted for 86.867% of the total original 
variance and most variables had high positive loadings, 
suggesting the size variation in the samples. The second 
principal component (PC2) accounted for 5.284% of the 
total variance and is strongly correlated with TBL, with 
factor loadings>0.9. Along the third principal 
component factor (PC3), which accounts for 3.863% of 
the total variance, the main morphological variable 
contributing to this association is PORCL, with factor 
loadings>0.9. 

Because the eigenvalue of the first principal 
components was 12.161, and the third was much lower 
(0.541), we used PCI and PC2 to make the PCA 
scatterplot (Figure 2). 


Table 1 Cranial variable measurements of Hylopetes (mean+SD)/range. Variable codes explained in Figure 1 





























; H. alboniger H. electilis H. nigripes H. platyurus H. spadiceus H. phayrei 
Уашаыз п=15 : n-22 kb 1 H Felt 0 ib n=19 

CRANL 48.52+3.02 38.56+0.99 51.50+2.70 31.90+0.88 35.58+0.80 40.65+2.56 
42.79—52.9 36.41— 39.96 45.59—55.28 30.87-33.26 34.58-36.87 37.51-46.61 

BCASEL 35.24+3.59 30.32+0.93 41.08+1.37 25.65+0.84 29.12+0.58 32.44+2.13 
30.17-41.4 28.47-31.61 38.06—42.94 24.20-27.00 28.35-29.72 30.04-37.61 

CRANW 21.98+0.62 18.21+0.44 22.58+0.52 16.38+0.39 17.67+0.29 18.84+1.12 
20.97-22.98 17.44-18.93 21.93-23.46 15.81-17.02 17.33-18.35 17.13-20.91 

BPORW 30.56+1.83 23.50+0.77 32.10+0.91 19.55+0.53 21.79+0.66 24.91+1.80 
26.84—33.3 22.23-24.74 30.14-33.06 18.77-20.25 20.98-23.04 21.97-28.34 

PORCL 11.51+0.79 11.20+0.65 11.71+0.77 9.78+0.52 10.20+0.57 10.10+0.64 
10.07-12.89 9.84—12.41 10.06-13.10 9.23—10.63 9.38-11.26 9.03-11.50 

PGA 19.63+1.37 15.42+0.60 20.79+0.96 13.66+0.44 15.40+0.64 16.65+1.20 
17.11-21.34 14.35-16.17 18.40—22.00 13.12-14.58 14.20—16.64 15.43-19.07 

NAL 14.35+1.36 10.66+0.72 16.17+0.85 9.06+0.41 10.15+0.48 11.67+1.13 
11.90-16.20 9.27-12.22 14.14-17.42 8.18—9.53 9.37—10.93 9.97—13.94 

TBL 8.79+0.52 8.55+0.51 9.30+0.46 6.64+0.30 7.40+0.17 9.43+0.47 
8.20—10.16 7.52—9.25 8.63—9.96 6.31—7.21 7.12—7.59 8.45—10.14 

DSL 10.21+1.06 8.15+0.40 10.46+0.77 6.42+0.38 7.33+0.40 8.79+0.77 
8.12-11.64 7.42—8.95 8.53-11.81 5.97—6.98 6.39—7.80 7.53-10.48 

MTRL 10.41+0.73 7.98+0.24 12.40+0.50 5.92+0.23 7.47+0.32 8.49+0.51 

8.00-11.06 7.45-8.42 11.77-13.18 5.62—6.28 7.01-8.19 7.73—9.82 

MTRW 12.86+0.54 9.96+0.28 12.96+0.29 7.74+0.27 8.66+0.40 10.33+0.56 
11.89-13.97 9.35-10.45 12.46-13.44 7.32—8.03 8.27—9.55 9.42-11.60 

LMDL 28.69+2.13 20.63+0.78 32.28+2.07 18.51+0.77 21.25+0.79 23.53+2.84 
25.65-33.57 19.42-21.97 28.32-35.33 17.48-19.73 20.26—22.60 18.46-28.00 

LMDH 18.42+1.56 14.22+0.62 21.05+1.20 11.38+0.64 13.02+0.64 14.82+1.09 
14.41-20.22 13.36-15.42 19.06-22.80 10.38-12.34 12.01-14.01 13.61-17.55 

LMTL 10.18+0.41 7.55+0.23 11.77+0.66 5.49+0.22 6.70+0.20 8.26+0.44 

9.50—10.87 7.19-8.04 11.04-13.56 5.18—5.82 6.46—7.13 7.72—9.36 
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Figure 2 indicates PCI-PC2 scatter plots of 
principal component analysis. Inspection of Figure 2, 
samples of H. electilis, H. platyurus, H. spadiceus, and H. 
phayrei are distinguished from each other to form their 
distinct groups, showing different morphological 
characteristics of their skulls. However, results showed 
that samples of H. alboniger and H. nigripes are quite 
close to one another with considerable overlap, 
suggesting that H. alboniger shares similar measurable 
skull characteristics with H. nigripes. An interesting 
result was that the three samples named Н. electilis from 
Chiang Mai, Thailand (labeled as “+” in Figure 2) are 
separated from those from Namfong, Hainan Island, but 
scatter in the samples of H. phayrei, hinting that the 
samples of H. electilis from Chiang Mai, Thailand should 
be more similar with those of H. phayrei in measurable 
skull characteristics. 





DISCUSSION 


Based on the specimens from Namfong, Hainan 
Island, China, Allen (1925) first described a valid species 
named Pteromys (Petinomys) electilis, but later on 
review classified it as Pteromys electilis (Allen, 1940). 
Later studies by Ellerman (1940) and Elleuman & 
Morrison-Scott (1950) continued to use Petinomys 
electilis as the scientific name for this species, and 
emphasized the species distribution on Hainan Island, 
China. However, based on their understanding of its two 
septa in each bulla, Corbet & Hill (1992) listed the group 
in Hylopetes. While Corbet & Hill (1992) narrated 
significant differences on pelage characteristics between 
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Н. electilis and Н. phayrei, they still opted to list it as 
belonging to H. phayrei, while placing H. p. electilis as 
its subspecies. More succinctly put, the validity of the 
Hainan flying squirrel, H.  electilis, has been 
controversial for quite a long time (Corbet & Hill, 1991, 
1992; Ellerman, 1940; Hoffmann et al, 1993; Nowak, 
1999; Thorington & Hoffmann, 2005; Wang, 2003; 
Zhang et al., 1997). 

The validity of various measurable cranial variables 
in distinguishing between species and species groups in 
Hylopetes has been critical in several studies. For 
example, Rasmussen & Thorington (2008) discussed the 
distributions for H. platyurus, H. spadiceus, and H. 
Lepidus using similar methods. In the present study, we 
performed statistical analyses on 14 measurable cranial 
variables from 22 samples of H. electilis, 19 of which 
came from Namfong, Hainan Island, China, in order to 
compare them with the samples of other five other 
validated Hylopetes species. The results showed 
significant difference between H. electilis and the other 
five Hylopetes species (Table 3), indicating it could be 
distinguished from other five species by using 
measurable skull variables. Similarly, though Figure 2 
indicates that samples of Н. electilis from Namfong, 
Hainan Island could be clearly distinguished from those 
of other five Hylopetes species, we found that they were 
closest to samples of H. phayrei, as opposed to the other 
four Hylopetes species. This result not only sustains the 
idea that H. electilis on Hainan Island should be a valid 
species, but it also may explain the reasons why there 
were so many systematic disagreements between H. 
electilis and H. phayrei. 


Table 2 Tests of Equality of Group Means (by Wilks' Lambda) for Sexual Dimorphism 





Wilks' Lambda F dfl df2 P 

CRANL 0.998 0.193 1 81 0.662 
BCASEL 0.996 0.332 1 81 0.566 
CRANW 1.000 0.022 1 81 0.884 
BPORW 0.999 0.074 Ji 81 0.787 
PORCL 0.995 0.443 1 81 0.507 
PGA 0.999 0.090 1 81 0.765 
NAL 1.000 0.038 1 81 0.846 
TBL 1.000 0.000 1 81 0.991 
DSL 0.996 0.308 1 81 0.580 
MTRL 0.998 0.154 1 81 0.696 
MTRW 1.000 0.008 1 81 0.930 
LMDL 0.998 0.148 dt 81 0.701 
LMDH 0.992 0.645 1 81 0.424 
LMTL 0.998 0.159 1 81 0.691 


Variable codes are explained in Figure 1. Significant difference level: ** 0< P<0.001, * 0.001< P<0.05. 
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Table 3 Multiple Comparisons (by Least-Significant Difference) of all 14 cranial variables between the six Hylopetes species 














Н. alboniger ~ Н. alboniger ~ Н. alboniger ~ H. alboniger — Н. alboniger ~ 
H. electilis H. nigripes H. platyurus H. spadiceus H. phayrei 
Mean Mean Mean Mean Mean 
difference P difference P difference P difference P difference P 

CRANL 0.097 0.000** —0.029 0.001** 0.180 0.000** 0.131 0.000** 0.075 0.000** 
BCASEL 0.063 0.000** —0.068 0.000** 0.136 0.000** 0.081 0.000** 0.034 0.000** 
CRANW 0.081 0.000** —0.011 0.072 0.129 0.000** 0.094 0.000** 0.067 0.000** 
BPORW 0.112 0.000** —0.022 0.006* 0.193 0.000** 0.146 0.000** 0.088 0.000** 
PORCL 0.014 0.120 —0.006 0.529 0.071 0.000** 0.054 0.000** 0.058 0.000** 
PGA 0.103 0.000** —0.024 0.008* 0.156 0.000** 0.106 0.000** 0.071 0.000** 
NAL 0.126 0.000** —0.053 0.000** 0.197 0.000** 0.149 0.000** 0.089 0.000** 
TBL 0.010 0.165 —0.026 0.003* 0.120 0.000** 0.072 0.000** —0.031 0.000** 
DSL 0.096 0.000** —0.010 0.418 0.200 0.000** 0.142 0.000** 0.065 0.000** 
MTRL 0.113 0.000** —0.077 0.000** 0.242 0.000** 0.142 0.000** 0.086 0.000** 
MTRW 0.109 0.000** —0.003 0.591 0.221 0.000** 0.172 0.000** 0.095 0.000** 
LMDL 0.142 0.000** —0.051 0.000** 0.187 0.000** 0.130 0.000** 0.087 0.000** 
LMDH 0.111 0.000** —0.059 0.000** 0.208 0.000** 0.150 0.000** 0.094 0.000** 
LMTL 0.130 0.000** —0.062 0.000** 0.270 0.000** 0.182 0.000** 0.092 0.000** 
CRANL —0.126 0.000** 0.083 0.000** 0.034 0.000** 0.021 0.001** 0.209 0.000** 
BCASEL —0.131 0.000** 0.073 0.000** 0.018 0.062 0.028 0.001** 0.205 0.000** 
CRANW —0.092 0.000** 0.047 0.000** 0.012 0.049* 0.014 0.005* 0.140 0.000** 
BPORW —0.135 0.000** 0.081 0.000** 0.034 0.000** —0.023 0.001** 0.216 0.000** 
PORCL —0.021 0.033* 0.056 0.000** 0.039 0.000** 0.044 0.000** 0.077 0.000** 
PGA —0.128 0.000** 0.052 0.000** 0.002 0.800 —0.032 0.000** 0.180 0.000** 
NAL —0.179 0.000** 0.070 0.000** 0.022 0.082 —0.037 0.001** 0.250 0.000** 
TBL —0.037 0.000** 0.109 0.000** 0.061 0.000** —0.042 0.000** 0.146 0.000** 
DSL —0.107 0.000** 0.103 0.000** 0.045 0.001** —0.031 0.004* 0.210 0.000** 
MTRL —0.191 0.000** 0.128 0.000** 0.028 0.002* 0.027 0.000** 0.320 0.000** 
MTRW —0.113 0.000** 0.111 0.000** 0.062 0.000** 0.014 0.013* 0.225 0.000** 
LMDL —0.194 0.000** 0.045 0.001** —0.011 0.362 0.055 0.000** 0.239 0.000** 
LMDH -0.171 0.000** 0.096 0.000** 0.038 0.001** 0.017 0.055 0.267 0.000** 
LMTL —0.192 0.000** 0.139 0.000** 0.051 0.000** 0.038 0.000** 0.332 0.000** 
CRANL 0.160 0.000** 0.104 0.000** —0.049 0.000** —0.105 0.000** —0.056 0.000** 
BCASEL 0.150 0.000** 0.103 0.000** —0.055 0.000** —0.102 0.000** —0.047 0.000** 
CRANW 0.105 0.000** 0.078 0.000** —0.035 0.000** —0.062 0.000** —0.027 0.000** 
BPORW 0.169 0.000** 0.111 0.000** —0.047 0.000** —0.105 0.000** —0.058 0.000** 
PORCL 0.060 0.000** 0.065 0.000** —0.017 0.162 —0.012 0.239 0.004 0.641 

PGA 0.130 0.000** 0.096 0.000** —0.050 0.000** —0.084 0.000** —0.034 0.000** 
NAL 0.202 0.000** 0.142 0.000** —0.048 0.002* —0.108 0.000** —0.060 0.000** 
TBL 0.098 0.000** —0.004 0.552 —0.048 0.000** —0.151 0.000** —0.103 0.000** 
DSL 0.152 0.000** 0.075 0.000** —0.058 0.000** —0.134 0.000** —0.076 0.000** 
MTRL 0.220 0.000** 0.164 0.000** —0.100 0.000** —0.155 0.000** —0.055 0.000** 
MTRW 0.176 0.000** 0.099 0.000** —0.049 0.000** —0.126 0.000** —0.077 0.000** 
LMDL 0.182 0.000** 0.138 0.000** —0.057 0.000** —0.100 0.000** —0.043 0.001** 
LMDH 0.209 0.000** 0.154 0.000** —0.058 0.000** —0.113 0.000** —0.055 0.000** 


LMTL 0.244 0.000** 0.154 0.000** —0.088 0.000** —0.177 0.000** —0.089 0.000** 
Variable codes are explained in Figure 1. **: 0< P « 0.001, *: 0.001« P « 0.05. 
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Table 4 Factor loadings of the principal component 








analysis 
Variables РС1 РС2 PC3 
CRANL 0.869 0.389 0.288 
BCASEL 0.851 0.403 0.223 
CRANW 0.875 0.262 0.346 
BPORW 0.872 0.345 0.313 
PORCL 0.295 0.112 0.945 
PGA 0.910 0.307 0.214 
NAL 0.912 0.284 0.214 
TBL 0.373 0.915 0.107 
DSL 0.775 0.523 0.214 
MTRL 0.844 0.363 0.304 
MTRW 0.801 0.416 0.367 
LMDL 0.935 0.222 0.174 
LMDH 0.886 0.326 0.279 
LMTL 0.819 0.417 0.321 
Eigenvalues 12.161 0.740 0.541 
Variance explained (%) 86.867 5.284 3.863 


Variable codes are explained in Figure 1. 


Both a variety of previous studies as all as the 
present study have established that the Hainan flying 
squirrel, H. electilis, is restricted to Hainan Island, China 
(Ellerman, 1940; Ellerman & Morrison-Scott, 1950), and 
we thought that the channel isolation may have 
accelerated its absolute evolution. Based on molecular 
data, Jorgensen & Demarais (1999) and Oshida et al 
(2000) discussed the close genetic relationship and the 
recent divergence between H. electilis and H. phayrei. 
Their results suggest that these populations rapidly 
extended the ranges to their present distributions in a 
relatively short time (Jorgensen & Demarais, 1999; 
Oshida et al, 2000). Based on our current findings and 
the previous data, we recommend further correlative 
ecology studies on H. electilis and H. phayrei should be 
performed in order to gather more data that would help 
explain their adaptation process, and in doing so explain 
their phenotypic differentiation. 

Additionally, three samples of Н. electilis (labeled 
as “+” in Figure 2) from Chiang Mai, Thailand, were 
found to be not only clearly different from those of H. 
electilis from Namfong, Hainan Island, China, but they 
also scatter within samples of H. phayrei to form a close 
group (Figure 2). Based on the results of pelage 
comparison between H. electilis and H. phayrei, while 
Corbet & Hill (1992) listed the former as H. p. electilis 
while also emphasizing some significant differences in 
pelage characteristics between the two (Corbet & Hill, 
1992). Our results also indicate that samples of H. 
electilis from Hainan Island could differ from those of H. 
phayrei samples from mainland in terms of skull 
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measurable characteristics (Figure 2), and likewise 
support their earlier results that H. electilis is distributed 
across Hainan Island, China. (Ellerman, 1940; Ellerman 
& Morrision-Scott, 1950). 

The ultimate results from our statistical analysis of 
the skull characteristics sustains the argument that Н. 
electilis is neither a synonym nor subspecies of H. 
phayrei, but a valid species in its own right. As for H. 
platyurus, although it was accepted as the synonyms of 
H. lepidus (Corbet & Hill, 1992), Figure 2 indicates that 
specimens of Н. platyurus could be completely separated 
from five other Hylopetes species, similarly illustrating 
that it likely must also belong to a valid Hylopetes 
species. Unfortunately, without specimens of H. lepidus 
to serve as a comparison, the data collected in the present 
paper is suggestive but not conclusive, therefore we 
cannot definitively state whether H. platyurus is a 
synonym of Н. lepidus or not. Further studies with more 
specimens may help further settle the ongoing debates 
regarding the taxonomy. 

For a final note regarding geography, it is worth 
noting that Hylopetes are distributed in different 
zoogeographic subregions of the Oriental region, each with 
significant geographical variations (Corbet & Hill, 1992; 
Lekagul & Mcneely, 1988; Zhang et al, 1997). Accordingly, 
the samples of the six Hylopetes species used in this study 
can be morphologically distinguished as different, distinct 
groups (Figure 2). H. alboniger and H. phayrei belong to 
the Indochinese subregion and overlap extensively in 
Thailand, Indochina, Burma, and southern China (Corbet 
& Hill 1992; Hoffmann et al, 1993; Thorington & 
Hoffmann, 2005), where the habitat is tropical, subtropical, 
or evergreen forest, with low elevation (500-1500 m) and 
within reach of southwest and south monsoons. 
Additionally, the population of Н. alboniger occurs in the 
eastern extreme of the Himalayas and occupies temperate 
coniferous forest at mid elevations (1 500-3 000 m) or 
deciduous and subtropical forests at low elevations, while 
H. nigripes is endemic to the Philippines and occurs in 
subtropical or tropical dry forests (Nowak, 1999). Similar 
morphological characteristics (Figure 2) between H. 
alboniger and H. nigripes may potentially be adaptations to 
similar living habitats, though again, further studies are 
needed to explore this possibility. 

South China itself is located at the crossroads of 
southeastern Asia, and has been a thoroughfare for 
animal dispersal from mainland Asia southward into the 
Indo-Malayan region. In southeastern Asia, many species 
have migrated south from Burma and southern China 
along the forested mountains of the Thai-Burma border 
into Malaysia, Sumatra, Java, and Borneo (Lekagul & 
Meneely, 1988). H. nigripes, H. lepidus, and Н. 
spadiceus inhabits the Sundaic sub region, with 
subtropical or tropical dry forests. The greatest 
distinction observed between H. platyurus and Н. 
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spadiceus and other Hylopetes is consistent with 
discussions about Н. platyurus, Н. spadiceus, and Н. 
lepidus by Rasmussen & Thorington (2008). However, 
because the habitats of these Sundaic Hylopetes are not 
well known and the ranges of these flying squirrels are 
poorly recorded, nearly every new collection will reveal 
new distribution limits, prompting further research. 
Based on that reality and the results of the present study, 
we suggest that future molecular studies on these species’ 
phylogenies may give us more evidence as to the volition 
of their species status, and the subsequent estimation of 
divergence times will help us to understand the 
biogeography of Hylopetes, and perhaps give some 
further insight into larger shifts within the region. 
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Specimens examined in this study 
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Yunnan, China: IOZ 5507 9, 5508 ĝ, 15046 9; KIZ 00009 9, 
76311 d, 84756, 89011, 90004, 206746, 206745, 
640222 &, 640227 3; AMNH 114884 3, 114885 3, 
114886 9; 
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Abstract: Sinocyclocheilus is a cave-dwelling cyprinid genus endemic to southwest China. Several species possess a conspicuous 
horn on their head, which has been suggested as a constructive troglomorphic trait but lacks substantial evidence. We used non- 
invasive, high spatial resolution synchrotron X-ray microtomography to investigate the three-dimensional (3D) morphology of the 
horn of Sinocyclocheilus hyalinus, one of eight such troglobiotic species. 3D renderings demonstrated the osteological components, 
which were comprised of a rear wall comprised of the supraoccipital bone, a remaining frontal wall with numerous fenestrae, and the 
bottom continuous with the parietal and epiotic. A horn cavity occurred within the horn. The fenestrae in the frontal wall were 
continuous in the horn cavity and showed elaborate channeling, and were, connected to the cranial cavity by soft tissue. We 
tentatively called this configuration the “otocornual connection” due to its anatomic and putative functional similarity to the 
otolateralic connection in clupeids and loricariids, which provide an indirect pathway to enhance perception of underwater sound 
signals. This study provides a functional morphology context for further histological and physiological investigations of such horn 
structures in Sinocyclocheilus cavefish, and we suggest that the horn might enhance acoustic perception to compensate for visual loss 


in subterranean life, which warrants future physiological examination as lab-reared S. hyalinus become available. 


Keywords: Sinocyclocheilus hyalinus, Cavefish; Horn; Troglomorphism; Synchrotron X-ray microtomography 


Abbreviations: C: cerebellum; Cc: cranial cavity; EPO: epiotic; F: 
frontal; Fo: foramen occipital; Foem: foramen occipital ellipticum magnum; 
FW: front wall of horn; He: horn cavity; LI: inferior lobe; MO: medulla 
oblongata; M: pd, muscle protractor dorsalis; OP: opercula; ot: optic tectum; 
P: parital; Pc: pharyngeal cavity; PG: pharyngeal bone; PT: pharyngeal teeth; 
PTE: pteroticSO: supraoccipital; TP: tripus; VC: valvula cerebella. 


Cyprinid genus Sinocyclocheilus is endemic to 
China with more than 50 valid species (Zhao & Zhang, 
2009). АП species exhibit cave-dwelling behavior to 
different degrees. To survive in such extreme subterranean 
conditions, cavefish have evolved a suite of morphological 
and physiological characteristics (troglomorphic traits) 
in both regressive and constructive directions, such as 
loss of eyes and pigmentation and enhancement of 
mechanosensory and/or chemosensory perception (Jeffery, 
2008, Poulson, 2010; Wilkens, 2003). A sub-group of 
Sinocyclocheilus cavefish possess a conspicuous horn in 
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their head, which is unique among known cavefish 
worldwide. The horn has been suggested as a 
troglomorphic trait but without substantial evidence to 
support this (Romero et al, 2009). In addition, little is 
known in regards to its anatomy or function, partially 
due to the limited availability of such rare specimens. 
Sinocyclocheilus hyalinus Chen, Yang & Zhu (1994), 
“Той тїр jin xian ba (translucent golden-line barbel)", 
is one such horn-containing, anophthalmic troglobiotic 
species (Chen et al, 1994). It was first reported from two 
specimens of 87 mm and 83.5 mm standard length, 
respectively. However, no further research has been 
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reported in the past two decades. X-ray tomography has 
been applied in large and medium-sized fish osteological 
studies. The advent of synchrotron X-ray microtomography 
expands this application into fine structures and soft 
tissue morphology without any staining in millimeter- 
centimeter-sized specimens (Pasco-Viel et al, 2010; 
Boistel et al, 2011). Synchrotron X-rays have the 
advantage of high coherence suited to phase contrast 
imaging, an intense photon flux, and a monochromatic 
and parallel beam with micron/sub-micron spatial 
resolution. 

In the present study, we employed synchrotron X- 
ray microtomography with 1.85-micron resolution to 
investigate the horn structure of one specimen of S. 
hyalinus. The resulting three-dimensional (3D) rendering 
and virtual section revealed the functional morphology of 
the horn. The resemblance of the horn configuration in 
the head of S. hyalinus and the “otolateralic connection” 
(such as recessus lateralis) in clupeids (Di Dario, 2004; 
Grande, 1985; Wilson et al, 2009) and loricariid catfish 
(Bleckmann et al, 1991; Schaefer & Aquino, 2000) 
suggest the possibility that the horn in Sinocyclocheilus 
fish might involve the enhancement of acoustic (and/or 
mechanosensory) perception. 


MATERIALS AND METHODS 


Specimen 

One specimen of S. hyalinus measuring 41 mm 
standard length (Figure 1) was collected in February, 
2009 from the Alu Ancient Cave (Luxi County, Yunnan 
Province), and was preserved in 10% formalin solution 
for about two years before the experiments. For the 
experiment, the specimen was switched to 70% ethanol 
solution and further dehydrated by gradient ethanol. A 
little shrinkage occurred after dehydration as indicated in 
the radiography (Figure 1). The horn was situated in the 
dorsal border between the head and trunk and displayed a 
forward protrusion. 


A 





Figure 1 Radiography (A) and photography (B) of the 
examined specimen of S. hyalinus 


Scale bar: 5 mm in A. One grid: 1 mm in B. 
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Synchrotron X-ray microtomography 

The experiments were performed at the X-ray 
imaging beamline (BLI3W1) at the Shanghai 
Synchrotron Radiation Facility (SSRF) in January and 
February, 2012. The S. hyalinus specimen was imaged 
by monochromatic synchrotron X-ray at 13.0 kV. The X- 
rays penetrating through the specimen were recorded as 
two-dimensional (2D) radiograph images (projection) by 
a charge coupled device (CCD) X-ray camera. The 
sample-to-detector distance was optimized for obtaining 
high in-line phase contrast. For synchrotron radiation X- 
ray microtomography (SR-wCT) scanning, a PCO.2000 
camera coupled with an optical magnification system 
(PCO AG, Kelheim, Germany) was used. In total, 720 
projections, with 2s exposure time for each, were 
recorded during specimen rotation over 180°. Flat-field 
images and dark-field images were also collected during 
each acquisition procedure to correct electronic noise and 
variations in X-ray source brightness. The data were 
reconstructed into corresponding slices using the filtered 
back-projection algorithm (X-TRACT software, CSIRO). 
The 3D renderings were created from the stack of slices, 
and were manipulated and analyzed in VG Studio Max 
(v2.1) software. 

During scanning, the specimen was held in a plastic 
tube vertically mounted on a sample stage. The 
synchrotron X-ray beam size was 45 mm (width) by 5 mm 
(height), so a 5 mm high segment including the horn was 
scanned using the 2X optical objective within the 
PCO.2000 camera. The local part containing the horn was 
then scanned using the 4X optical objective covering a 
height of 3.78 mm. The resulting spatial resolutions were 
3.7 micron/pixel and 1.85 micron/pixel, respectively. 

Additionally, a 2D radiography of the whole 
specimen (Figure 1) was stitched by two consecutive 
radiographies recorded using a VHRI:1 X-ray digital 
camera with direct coupled (micro) fiber-optic with 9 
um/pixel resolution (Photonics Science Ltd, 
Robertsbridge, UK). The specimen was conventionally 
positioned in soft plastic foam on a sample stage, which 
vertically moved approximately 5 mm to allow the two 
consecutive radiographies. 


RESULTS 


Reconstructed 2D slice and 3D rendering 

The original reconstructed slice from the SR-uCT 
was an 8-bit grayscale image on the coronal plane. The 
mineralized tissue (pharyngeal teeth, vertebra, other 
bones), soft tissue (brain, musculature), and ambient air 
(cranial cavity, pharyngeal cavity) were easily distinguished 
from one another according to their different gray values 
(Figure 2A, C). The mineralized tissues with high grey 
values were white in the reconstructed slices, whereas air 
with the lowest value was dark. The 3D rendering was 
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created in VG Studio software using the stack of slices, 
from which sagittal and horizontal virtual sections could 
be further created (Figure 2C). By increasing grey pixel 
value, the pixel below or equal to the set threshold was 
interpreted as black. Such manipulation allowed for the 





virtual removal of the soft tissue (“virtual corrosion”), 
and demonstrated the osteological components (Figure 
2B, D-I). Areas with the same grey value could be 
delimitated and labeled in designated colors, thus the 
potential cavity inside the specimen could be labeled. 


crests — 


SO cre 


Figure 2 2D reconstructed slices and overall 3D rendering of the horn 


Each reconstructed coronal and sagittal slice is presented (A, C), where bones, soft tissue, and air-filled cavities are distinguished from each other by different 


grey values. Overall 3D rendering (B) was created from the stack of original reconstructed coronal slices. After virtual removal of musculature and other soft 


tissue, the osteological components are displayed in lateral, oblique, and rear view (0-Е). A—F are from SR-uCT scanning at 3.7-miron resolution. G-I are 


from scanning at 1.85-micron resolution. The top part of the horn remains in G, but is removed in H and I to display the bony bottom of the horn. The horn 


cavity is labeled in pale green. H and I differ in whether the remaining horn cavity is labeled or not. A star marks the fenestrae on the frontal wall. Scale bar: 1 


mm in A; 1.5 mm in B, D-F; 0.85 mm in C, G-I. 
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General 3D morphology of the horn 

To observe the overall 3D morphology of the horn, 
a 3D rendering of a 5-mm segment containing the otic 
region and the horn was created (Figure 2B). The horn 
was situated at the dorsal otic region of the skull and bent 
forward as a protrusion with a flattened ending. After 
virtual removal of musculatures and other soft tissues, 
the osteological components of the horn were presented, 
which were divided into rear wall, remaining frontal wall, 
and bottom. The supraoccipital (SO) bone provided the 
rear wall (Figure 2F), at which the hypertrophy muscle 
protractor dorsalis (M.pd) was attached (Figure 2C). The 
supraoccipital crest and two additional crests on both 


sides formed along the suture between the SO and 
dorsal-posterior parts of the remaining frontal wall (fW) 
(Figure 2F). The fW, bearing numerous fenestrae, arose 
from the parietal and epiotic (Figure 2D, E). The bottom 
was continuous with the parietal and the epiotic, but we 
were unable to identify whether the parietal or the epiotic 
formed the bottom. The SO and the fW merged at the 
distal ending and separated at the basal part of the horn. 
Thus, a cavity was formed within the horn, and the 
channeling of the fenestrae in fW connected to it was 
more elaborate. The horn cavity was labeled in pale 
green in individual panels (Figures 2 and 3), and further 
demonstrated the following. 





Figure 3 Connecting pathway “frontal wall bearing elaborate fenestrae horn cavity with soft tissue connection cranium" revealed 


in representative horizontal (B—D), coronal (F-H), sagittal (J—L) sections 

In the two upper horizontal sections, the horn cavity with trabeculae-like structure and the fenestrae with elaborating channeling are 
most evident (B, C), while in the lowest horizontal representative section the cranium appears and the bony bottom with one foramen 
is also shown (D). In the coronal sections, the cranial cavity containing part of cerebellum is evident, and elaborate channeling of 
fenestrae is also visible (F, G); the horn cavity almost disappears in the posterior representative section (H). The fenestrae and 
foramen are indicated with a star. In sagittal sections, the soft tissue, which arose from the dorsal medulla oblongata and connected 
the horn and cranial cavities, are most evident, and the foramen position through which the connection penetrated are indicated by 
arrows (J-K). Skin tissue investing the horn is also displayed (Figure 3B-D, F, О, J). In each row, the position of the representative 
sections are indicated in the left most panels, where the horn cavity labelled in pale green is visible through the fenestrae in fW or the 
foramen in the bony bottom (A, E, I). All images are from SR-uCT at 1.85-miron resolution. Scale bar: 1.1 mm in A,E, I; 0.55 mm in 
B-D, J-L; 0.4 mm in F-H. 
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Configuration of “otocornual connection" 

To better demonstrate details, a 3D rendering from 
an additional SR-uCT scan at 1.85-micron resolution 
was created. In the oblique view, the horn cavity from 
the cranial cavity was evident after removal of the top 
part of the horn (Figure 2G-D, especially through 
comparison between the images with and without 
labeling of the remaining horn cavity (Figure 2H-I). 
Several foramens in the bony bottom were evident from 
the ventral view, through which the horn was connected 
to the cranial cavity (Figure 3I). 

Furthermore, sagittal, coronal, and horizontal virtual 
sections were performed. A configuration, tentatively 
called the *otocornual connection", was resolved from 
these sections. In general, the fenestrae in fW showed 
elaborate channeling and were continuous with the horn 
cavity , which in turn was connected with the cranial 
cavity by soft tissue. The horn cavity and fenestrae with 
elaborate channeling were most evident in the 
representative horizontal sections (Figure 3B-D), but 
were also easily identified in the sagittal and coronal 
sections where the fenestrae were labeled with a star 
(Figure ЗЕ, 3G, 3J-L). The soft tissues forming the 
connection between the horn and cranial cavity were 


0.65 mm 





evident in the representative sagittal sections, which 
arose from the dorsal medulla oblongata and penetrated 
into the horn cavity through the bony bottom (indicated 
by arrows in Figure 3J-L). Unfortunately, we were 
unable to further identify whether the soft tissue 
connecting the two cavities was neural or membrane 
tissue based on the grey value in the image. The 
trabeculae-like soft tissues (or connective tissue) were 
also visible in the horn cavity in the horizontal and 
sagittal sections (Figure 3B, C, K, L). The cerebellum, 
otic tectum and part of the medulla of the brain with 
certain decay were also shown in the sagittal sections 
(Figure 3J—L). 

To give an intuitive presentation, two modified 3D 
renderings were created with the horn cavity labeled in 
pale green, and the cerebellum and otic tectum labeled in 
khaki-grey as a reference for the cranium (Figure 4). In 
one 3D rendering, the encasing structures were set in 
semitransparent status to emphasis the horn cavity and 
brain tissue (Figure 4A). In the other rendering, the left 
half encasing structure was removed to exposure the 
horn and cranial cavity, where the connecting soft tissue 
was visible and appeared to be membrane rather than 
fiber tissue (Figure 4B). 


0.65 mm 


Figure 4 Two 3D renderings for the intuitive presentation of the *otocornual connection" configuration 


The horn cavity is labeled pale green, and the cerebellum and otic tectum khaki-grey. The encasing structure is semitransparent in A, 


whereas the left half encasing structure is removed to exposure the horn cavity and cranium in B. The connecting soft tissue is 


indicated by an arrow in B. 


The “otocornual connection" might provide a 
“frontal wall bearing elaborate fenestrae — horn cavity 
with soft tissue connection—cranium” pathway, with its 
potential function discussed in next section. 


DISCUSSION 


We employed a non-invasive SR-uCT approach to 
investigate a S. hyalinus specimen at resolution of 1.85- 
microns, and revealed the osteological components and 
“otocornual connection” configuration of the horn. The 
3D morphology results will help clarify previous reports 
in relation to anatomic structure, and provide evidence 
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about the possible function of the horn to enhance 
acoustic perception to compensate for visual loss in cave 
life. 

The Sinocyclocheilus horn is unique among all 
known cavefish species. Eight troglobiotic species, 
namely S. angularis, S. aquihornes, S. bicornutus, S. 
broadihornes, S.  furcodorsalis, S. hyalinus, S. 
rhinocerous, and S. tileihornes, possess these horny 
protrusions. Wang et al (1999) delineated the horn of S. 
angularis and S. bicornutus as a whole parietal 
protrusion without supraoccipital or other osteological 
components, which was concluded in subsequent 
literature that the horn consisted of the frontal and 
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parietal bones. Li et al. (2002) performed gross anatomy 
on S. rhinocerous and concluded that the horn was 
composed of three ossicles set perpendicular to one 
another. Our 3D morphology results showed that the 
supraoccipital bone and ossicles arose from the parietal 
and epiotic, and a bony bottom comprised Ше 
osteological base of the horn. These results indicated that 
the frontal did not build the horn structure, and disagreed 
with previous findings that the horn is a “whole parietal 
protrusion”. The explanation that the three ossicles are 
set perpendicular to one other in S. rhinocerous seems to 
be an incorrect interpretation of the three crests 
(supraoccipital crest and two additional crests on both 
sides formed by suture of the supraoccipital and posterior 
edge of the front wall) found in the rear wall in S. 
hyalinus. Recently, molecular phylogenetic analysis has 
indicated that all examined horn-possessing species are 
clustered within one monophyletic clade (Xiao et al, 
2005; Li et al, 2008), called the “angularis” clade (Zhao 
et al, 2009). Systematic investigation of the 3D 
morphology of horns and analysis of the conservation 
and variation among them would benefit our 
understanding of evolution within this “angularis” clade. 

The horn has been suggested as а possible 
troglomorphic characteristic for Sinocyclocheilus 
cavefish (Romero et al, 2009), but whether this structure 
is related to adapted subterranean life, the function of the 
horn, and how it contributes to sensory perception are 
still unknown. From the 3D morphology, it is tempting to 
speculate that the horn might be an adaptation to enhance 
acoustic perception. Underwater sound has two physical 
components, particle motion and sound pressure wave. 
All fish can detect particle motion through inertial 
movement between the otolith and sensory epithelia 
within the inner ear. Some fish can also detect sound 
pressure by physical connection structures transforming 
sound pressure into the inner ear (Popper, 1993, 2011). In 
otophysans, the Weberian ossicles connect the swim 
bladder and otic capsule. The swim bladder inside the 
abdominal cavity is set into motion by pressure in the 
sound field, and such vibrations are transformed by 
Weberian ossicles into particle motion perceived by the 
inner ear. In clupeids, the complex recessus lateralis 
structure in the otic region of the skull is an intracranial 
space into which the supraorbital,  infraorbital, 
preopercular and temporal sensory canals open (Grande, 
1985; Di Dario, 2004). A pair of gas-filled auditory 
bullae in clupeid heads directly contacts the recessus 
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lateralis at the lateral side and the modified utricle at the 
medial side. Proper configuration of the recessus 
lateralis, auditory bullae and modified utricle provides 
an indirect pathway for the pressure component of sound 
to be perceived by the inner ear of certain clupeid species. 
Recessus lateralis lesions or filling auditory bullae with 
ringer solution has been shown to damage the ability to 
detect 40 kHz sounds in Gulf menhaden (Wilson et al, 
2009). In loricariid, a similar structure is also found in 
the temporal region of the skull, where a section of the 
posotic laterosensory canal deep within the ventral layer 
of the compound pterotic is in direct contact with the 
horizontal canal of the inner ear (Bleckmann et al, 1991; 
Schaefer, 2000). These "pressure transducer" otophysical 
connections in otophysans and otolateralic connections 
in clupeids improve hearing capability (such as auditory 
sensitivity, frequency bandwidth, and sound source 
location). In the present 3D rendering of S. hyalinus horn, 
а “frontal wall bearing elaborate fenestrae horn cavity 
with soft tissue connection cranium" configuration was 
clearly revealed, which was tentatively named the 
"otocornual connection" due to its anatomical and 
putative functional similarity to  otophysical and 
otolateralic connections, which enhance sound pressure 
perception. Acoustic perception could be particularly 
important in the subterranean lightless environment of 
blind cavefish. Whether the horn functions to enhance 
acoustic perception warrants future physiological studies 
as lab-reared S. hyalinus become available. 

In conclusion, we demonstrated the presence of a 
cavity inside the horn of S. hyalinus with a connection to 
the cranium. The supraoccipital bone provides the rear 
wall of the horn, and the remaining frontal wall bears 
numerous fenestrae. Such a configuration implies a 
potential putative function as an indirect pathway to 
perceive sound pressure. This study also provides a 
functional morphological context for further anatomic 
and functional investigations of horn structures found in 
Chinese Sinocyclocheilus cavefish. 
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Effects of Zn^* and Си“ on loach ovaries and ova development 
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Abstract: This study compared the accumulation of Zn” and Cu” in the ovaries and ova of loaches under different concentrations of 
Zn?" (1.00, 2.50 and 5.00 mg/L respectively) and Cu?" (0.10, 0.25 and 0.50 mg/L respectively). The results showed that both Zn™ and 
Cu? accumulated in the ovaries, and that the relationship between accumulation and time was linear over 20 days of exposure. The 


accumulation of the metals in ovaries was closely related to the concentration of exposure in the solutions (P«0.05), and was 


obviously affected by the time and doses. However, the Cu?” concentration was significantly higher than Zn?' (P«0.05). The 


development level of ova in the ovaries also correlated with the concentration and exposure period in the Zn^* and Cu” solutions. 


Keywords: Misgurnus anguillicaudatus, Ovary; Ovum; Zn; Cu; Accumulation 


The problem of heavy metal pollution in aquatic 
environments has become an increasingly serious 
concern over the past decades in China (Wang et al, 
2008). The resulting pollution that accompanies 
industrialization has severely challenged the survival of 
aquatic animals due to strong toxicity and the 
bioaccumulation of heavy metals. Of these pollutants, 
metal ions accumulate not only in the skin, muscles, liver, 
and kidneys (Al-Weher S et al, 2008; Dutta T et al, 2001; 
Has S et al, 2008; Migliarini et al, 2005), but also in the 
gonads, which poses a serious threat to the reproductive 
success of many aquatic species and by extension the 
overall population of those animals (Abou El-Naga et al, 
2005; da Cruz et al, 2007; Tang et al, 2010). Accordingly, 
to maintain biodiversity and protect the breeding system 
of aquatic fauna, it is necessary to conduct studies on the 
accumulation of heavy metals in a fish's reproductive 
organs and the resulting stress on their ova to thereby 
establish correlations between metal accumulation and 
lon concentration. 

Zn?" and Cu” are commonly found in bodies of 
water polluted by heavy metals. Neither biodegrades 
easily and these metals are amplified via bioaccumulation 
as they pass through food chain. To explore how heavy 
metal accumulation strains the development of the ovaries 
and ova, in this study we used Loaches (Misgurnus 
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anguillicaudatus), a fish species known to tolerate 
pollution well (Gao et al, 2003; Wang et al, 2003) as a 
model. By doing so, we hope to find results that will 
prove useful in gauging the impact assessment of aquatic 
organisms, diagnostics of environmental pollution, and 
biodiversity conservation. 


MATERIALS AND METHODS 


Animal specimens and reagent materials 

Two year old loaches (n=220) were purchased from 
a local farmer's market for use as a testing species. The 
overall mean length of the fish was 13.2+3.5 cm and the 
mean body weight 3750.48 g. 

Concentrated solutions (1 000 mg/L) of ZnSO, (AR) 
and CuCl, (GR) were first prepared before being diluted 
into the corresponding concentrations as needed throughout 
the course of the experiment. HNO; (AR) and HCIO, 
(AR) were mixed at a ratio of 4:1 before use. 








Instruments and equipments 
The following equipment were used over the duration 
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of this study: an H-6800 Aerating Pump (China); an 
IRIS Intrepid ER/S-model ICP atomic emission 
spectrometer (Thermo Elemental Co., USA); a CKX41 
inverted fluorescent microscope (HQ2592x1944) 
(Olympus Co., Japan); an LSP far-infrared cooking 
stove (China); a KD (freeze) slicing machine (China); 
and a DGG-9070A drying oven with an electrical 
thermostatic wind drum (China). 


Experimental design 

The experiment was carried out indoors. During 
testing, the water quality parameters were maintained 
consistently: pH 6.3—6.5, DO 5.1—6.2 mg/L; temperature 
at 8-12°C (measured daily); average total water hardness 
at 2.67 mmol/L; and average alkalinity at 2.6 mmol/L. 

Polyethylene plastic aquariums (40 cmx30 cmx45 cm) 
were used as exposure devices, and 20 L of tap water, in 
which no Zn?* and Cu” was detected, was added to each 
tank and then aerated for more than 2 days prior to being 
used. The loaches accepted for testing were visually 
inspected to be disease-free, injury-free and active, with 
an overall mortality rate less than 5% during the testing 
period. The fish were acclimated to lab conditions for 5 
days before the experiment began. Over the course of 
testing, there was no water exchange in the plastic 
aquariums and the fish were not fed. An aerating pump 
was used to keep dissolved oxygen above 5 mg/L. 

In total, 220 similarly sized loaches were randomly 
placed into four testing groups. The control group was 
placed in water with no heavy metals present, while the 
other three were given different concentrations of Zn? 
and Си?" similar to actual ion concentrations in the local 
water bodies based on the Chinese Standard of Water 
Quality for Fisheries (GB11607-1989): Zn** 
concentrations of 1.00, 2.50 and 5.00 mg/L paired with 
Си?* concentrations of 0.10, 0.25 and 0.50 mg/L. The 
fish were held in each aquarium and three parallels for 
each concentration over the course of the 20-day 
experiment. Each individual was weighed, boiled, and 
dried prior to analysis, during which ovary slices from 


different groups and at different times (before testing, 5- 
day, 10-day and 20-day) were cut and observed. 


Statistical analysis 

The processing method of statistics for the testing 
results follows the method described by Nan Xu-Yang 
(Nan, 2009). All data was analyzed using SPSS 11.0 data 
processing system (SPSS Inc., Chicago, USA) for analysis 
of variance. The values for individual experiments were 
collected to calculate the mean value and the standard 
deviation/error to make comparisons. The statistical 
significance of the difference between means was 
determined using one-way ANOVA, with P<0.05 being 
Statistically significant. 


RESULTS 


Accumulating rule of Zn^* and Cu” in ovaries 

Before exposure to the metals, Zn?" and Cu”, all 
loaches were tested to ensure that no metals were present 
in the ovaries. By 10 days of exposure, the concentrations 
of Zn^* and Cu” in loach ovaries across all groups being 
kept in water with varying concentrations of heavy 
metals had increased. On day 20, the concentration of 
both ions showed an increasing trend, but between day 10 
and day 20 the concentrations of Zn^' and Си” in the 
ovaries seemed to decline significantly when compared to 
ion concentrations of the groups when they had been 
exposed for less than 10 days (P«0.05). 

The cumulative capacity of Zn^' and Си? in the ova 
was positively related to concentration of the metals in 
aqueous solution (P«0.05). The accumulation of Zn^ 
and Cu? in the ovaries presented as a function of 
exposure time and dose effects are shown in Table 1. In 
the control group, по Zn^' and Cu^* were detected. 

Table 2 shows the regression equations and 
correlation coefficients for ovaries, indicating a positive 
relationship between Zn" and Cu^ accumulation and 
exposure time. Relationships between accumulation of 
Zn^ and Cu^ in the ovaries and exposure time are 
shown in Figures 1 and 2. 


Table 1 Accumulation of Zn^* and Cu” in ovary at different time intervals (12200) 





























Time (d) Treatment Zn” / Cu in solution (mg/L) Zn” (ug/g; dry wt) Cu” (ug/g; dry wt) 

5 control — — — 
group 1 1.00 / 0.10 1.3+0.48 3.4+0.28 
group 2 2.50 / 0.25 2.5+0.4® 5.6+0.8® 
group 3 5.00 / 0.50 3.6+0,7®! 8.9+0.98 

10 control — — — 
group 1 1.00 / 0.10 2.3+0.7° 54.342.5! 
group 2 2.50 / 0.25 5.1+0.2° 81.5+4.7° 
group 3 5.00 / 0.50 26.441.7° 145.3345.6° 

20 control — — — 
group 1 1.00 / 0.10 4.3+0.3% 117.1+6.7° 
group 2 2.50 / 0.25 7.6+0.4° 137.2+11.0° 
group 3 5.00 / 0.50 39,143.28 329.3+18.4° 





Different superscripts show significant difference (P < 0.05) , while the same superscript shows no significance (P>0.05) between groups. 
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Table 2 Regression equations and coefficients of accumulation for Zn?* and Cu with exposure time in ovary 


Regression equation 


Correlation coefficient (R?) 





y=0.185x+0.7408 0.7834 
Y=9.3235x-30.158 0.9543 
y=0.3254x+1.2608 0.9355 
Y-9.3654x-15.307 0.8286 
y-2.2072x-2.1333 0.8648 

0.9716 


y=24.754x-87.721 







Treatment Heavy metal 
group 1 Zn 
Cu 
group 2 Zn 
Cu 
group 3 Zn 
Cu 
35 f —— 0.10 mgL" 
30 f —— 0.25 mgL" 


= 0.50 mg L” 


Cumulant of Zn 
(ug/gDW) 


5 10 20 
Time(d) 


Figure 1 Relationship of accumulation of zinc and time in ovary 


ES 0.10 mg: L” 


+ 0.25 mg: L`' 
-— 0.50 mg: L” 


Cumulant of Cu 
(ug/gDW) 
wo £ z Š Š 


5 10 20 
Time(d) 


Figure 2 Relationship of accumulation of copper and time in ovary 


Stress on ovary development from Zn” and Cu? 

The ovaries and ova in the controls, in which the 
metals were not detected, developed normally. Zn?" and 
Си?“ accumulated rapidly in the ovaries of all treatment 
groups exposed to mixtures of Zn?" and Cu”. The fish 
showed on the atrophy and cytoplasmic leakage in the 
ova, suggesting histological damage in the ovaries by 
day 10 of exposure. By day 20, the ova showed 
symptoms of severe degeneration, the mutual bonding of 
cells and severe atrophy (Figure 3). 


DISCUSSION 


As trace elements of an organism, Zn^' and Cu” are 
integral parts of cellular structure and enzyme 
compounds. Zn?" is one of the important growth factors 
and plays a vital role in growth, development and 
appetite, etc. for fishes and other aquatic organisms (Al- 
Weher S, 2008; Ebrahimi M, 2007). However, studies 
have shown that higher doses and longer exposure are 
not beneficial, and have toxic effects on aquatic 
organisms (Zhou et al, 2002). As for Cu”, it is also a 
type of nutritive trace element necessary for fish growth 
(Gao et al, 2003), and it performs important functions for 
the physiological activities in fishes. The physiological 
metabolism will be negatively influenced if either Хп?“ 
or Cu” is lacking, but overdosing brings adverse impacts, 
especially in the case of Cu’. In general, the 
harmfulness of Zn” is relatively small as compared with 
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Figure 3 Effects of joint attack of Zn?' and Cu?" 


on ovaries and ova 


a: Normal development of ovaries and ova in the control; b:5 d ovaries and 
ova, the arrow shows the cell degenerated in group one; c:10 d ovaries and 
ova, the arrow shows the cell dehydrated in group one; d:10 d ovaries and 
ova, the arrow shows the cell atrophic and dehydrated in group two; e:10 d 
ovaries and ova, the arrow shows the cell atrophic and shriveled in group 
three; f:20 d ovaries and ova, the arrow shows the cell atrophic and mutual 
bonding in group one; g:20 d ovaries and ova, the arrow shows the cell 
atrophic and mutual bonding in group two; h:20 d ovaries and ova, the 
arrow shows the cell atrophic and mutual bonding in group three. 


Cu^, which is highly toxic to aquatic organisms in 
stronger dosages (Wang et al, 2006). 

As mentioned above, the amount of Zn?" and Си“ 
accumulated in the ovaries was related to the 
concentration of Zn? and Cu? in the aquatic 
environment and the duration of exposure to the metals. 
The cumulative amount of heavy metal in the ovaries 
rises with the concentration of Zn?* and Cu” to a certain 
level alongside the exposure period, demonstrating the 
significant effect of time and doses. However, the 
accumulation rate tended to decline when exposure time 
exceeded 10 days. It could be that when organisms are 
over-exposed to heavy metals such as Zn?" and Cu”, the 
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metals activate the transcription of metallothionein genes 
in the organs, resulting in more genes being expressed. 
Consequently, these expressed genes combine the metals 
that have entered the cells with those new synthetic 
proteins (Allen, 1995). Alternatively the cells over- 
exposed to the heavy metals are in a state of “saturation” 
and thus can no longer absorb more heavy metals. As 
exposure continues, the level of heavy metal 
accumulation declined in the ova, which was more 
obvious in the cases of different concentrations of Zn^ 
(Figure 1). As for the accumulation of Си?“ under 
different concentrations, the higher the concentration of 
Cu”, the slower the accumulation tended to be. It 
seemed that the cumulated amount of Cu?' increased in 
the ova between 10 to 20 days, when the concentration 
was low, 1.е., 0.1 mg/L (Figure 2). 

Ultimately, these results indicate that the ability of 
the ovaries to detoxify decreased and the ova were 
damaged as a result of 20 day’s exposure to heavy metal 
at higher concentrations. Long term exposure to heavy 
metals will likely then further damage the overall 
function of ovaries (Figure 3), which is in agreement 
with previous studies regarding organs in aquatic 
organisms being damaged by heavy metals (Wang et al, 
2003; Zhang et al, 2009; Zhou et al, 2010). 
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Abstract: Pratylenchus ekrami from maize (Zea mays) roots in Shenyang and luffa (Luffa cylindrica) roots in Dalian, China, are 


described in this paper. Nematodes from the two areas were identified consistently, and were characterized by a heavy cephalic 


sclerotization, extending posteriorly up to two body annuli, stylet 11—13 шт long, elongating conoid tail, and becoming thinner from 


vulva. Males were not found. Pratylenchus ekrami is close to Pratylenchus vulnus, but the most critical characteristics between the 


two species were the number of lip annuli, stylet size, and shape of stylet knobs. This species is the first reported in China. 


Keywords: Pratylenchidae; Pratylenchus; new record species; Pratylenchus ekrami 


Root-lesion nematodes of the genus Pratylenchus 
(Bajaj & Bhatti, 1984) rank second only to root-knot and 
cyst nematodes in terms of worldwide economic impact 
on crops (Sasser & Freckman, 1987). This is due not 
only to their wide host range, but also their distribution 


in almost every cool, temperate and tropical environment. 


Pratylenchus species are obligate biotrophic, soil- 
inhabiting parasites and are found in all agricultural 
regions of the world (Hunt & Perry, 2007). So far, more 
than sixty species in the genus have been recorded 
(Castillo & Vovlas, 2006). Pratylenchus ekrami Bajaj & 
Bhatti, 1984, belonging to Tylenchida, Pratylenchidae 
has been recorded only in the soil and roots of Pyrus 
malus L. at a horticultural farm in India. Body of the 
genus is open C-shaped upon fixation and becomes 
thinner from the vulva; labial framework is sclerotized, 
extending posteriorly up to two body annuli; stylet is 
11-13 um long and the tail is conoid (Bajaj & Bhatti, 
1984). One species was found from maize (Zea mays) 
roots in Shenyang and luffa (Luffa cylindrica) roots in 
Dalian, China. It was identified as Р. ekrami after 
morphological studies, and is reported here for the first 
time. 


MATERIALS AND METHODS 


All specimens used in this study were obtained 
from the soil of maize roots in Shenyang and luffa roots 
in Dalian, Liaoning Province, which are the first 
reported occurences of this species in China. Soil 
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samples were dissociated by sucrose centrifugation 
(Chen et al, 2007). Specimens were fixed in FG (8 mL 
40% formaldehyde, 2 mL triethanolamine and 90 mL 
sterile water) for measurement and photographs. First, 
these specimens were placed on permanent slides. The 
shape of the females was then observed directly using 
optical and biological microscopes. The fine structure 
was observed by using an oil immersion lens. 
Photographs were taken and features were measured 
using motic metrology software. Good photos were 
selected to be drawn. 


RESULTS 


Pratylenchus ekrami Bajaj & Bhatti, 1984, a new 
record species in China (Figure 1). 


Measurements 

Female (n=10): L=422.0—490.0 (449.9425.20) um, 
W-13.7—20.2 (16.942.43) um, a—21.7—33.3 (26.4+4.46), 
b—5.2-6.7 (6.03+0.47), Б=3.4—5.2 (4.5+0.85), с 
=15.3-22.6 (18.88+2.34), с=1.6-22 (1.83+0.26), 
ү=79.9%—83% (81.58+1.07), stylet-11.2-13.1 
(12.23+0.74) um, Tail=18.7—28.3(23.4843.13) um. 
































Received: 05 March 2013; Accepted: 15 June 2013 

Foundation item: This work was supported by China Agriculture 
Research System (CARS-04) 

*Corresponding author, E-mail: duanyx6407@163.com 


Volume 34 Issues E4—5 


Description of a new Pratylenchus species from China (Tylenchida, Pratylenchidae) E141 


L: body length; W: greatest body width; a: body 


length/greatest body width; b: body length/distance from 
anterior end to junction of oesophagus and intestine; р": 
body length/distance from anterior end to oesophagus 
end; c: body length/tail length; с”: tail length/anus body 


width; V: distance from vulva to anterior endx100/body 
length. 


Description 

Female: Body opened C-shaped upon fixation. The 
body became slender from vuval region. Lateral fields 
with four lines extended to tail terminus. Labial region 
was continuous with the rest of the body, anteriorly 
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truncate, and had three annules. Labial framework was 
sclerotized, extending posteriorly up to two body annuli. 
Stylets were 11-13 um long with rounded, sloping 
knobs. Corpus was long, almost cylindrical. Median 
bulb was oval. Pharyngeal gland overlapped the 
intestine. Opening of dorsal was 3—4 рт from stylet 
base. Nerve ring was located at mid-isthmus. Excretory 
pore was 74—79 um from anterior extremity. Hemizonid 
was just anterior to excretory pore. Spermatheca was 
elongated oval and filled with sperm. Vulva was a 
transverse slit located at 77%—83% of body length. Tail 
was elongated conoid or blunt. 
Male: Not found. 
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Figure 1 Illustrations of Pratylenchus ekrami 
A: Female pharyngeal region; B: Tail of female; C: Female labial region; D: Female esophageal bulb region; E: Female vulvar region; F: Lateral field; G: Entire 


female. 


Kunming Institute of Zoology (CAS), China Zoological Society 


Volume 34 Issues E4—5 


E142 LIU, et al. 


DISCUSSION 


Pratylenchus ekrami was close to Pratylenchus vulnus, 
although they differed in shape of labial region, number of 
lip annuli, stylet size, shape of stylet knobs and median bulb, 
details of lateral lines and number of tail annuli. 
Pratylenchus ekrami is characterized by: labial region 
continuous with rest of body, anteriorly truncate, three 
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